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ABSTRACT 

We report direct evidence for asymmetry in the early phases of SN 1987 A via optical spectroscopy of five 
fields of its light echo system. The light echoes allow the first few hundred days of the explosion to be reob- 
served, with different position angles providing different viewing angles to the supernova. Light echo spec- 
troscopy therefore allows a direct spectroscopic comparison of light originating from different regions of the 
photosphere during the early phases of SN 1987 A. Gemini multi-object spectroscopy of the light echo fields 
shows fine-structure in the Ha line as a smooth function of position angle on the near-circular light echo rings. 
Ha profiles originating from the northern hemisphere of SN 1987 A show an excess in redshifted emission and 
a blue knee, while southern hemisphere profiles show an excess of blueshifted Ha emission and a red knee. 
This fine-structure is reminiscent of the "Bochum event" originally observed for SN 1987 A, but in an exag- 
gerated form. Maximum deviation from symmetry in the Ha line is observed at position angles 16° and 186°, 
consistent with the major-axis of the expanding elongated ejecta. The asymmetry signature observed in the 
Ha line smoothly diminishes as a function of viewing angle away from the poles of the elongated ejecta. We 
propose an asymmetric two-sided distribution of 56 Ni most dominant in the southern far quadrant of SN 1987A 
as the most probable explanation of the observed light echo spectra. This is evidence that the asymmetry of 
high-velocity 56 Ni in the first few hundred days after explosion is correlated to the geometry of the ejecta some 
25 years later. 

Subject headings: supernovae: individual(SN 1987A) — supernova: general — supernova remnants — line: 
profiles 



1. INTRODUCTION 

A light echo (LE) occurs when outburst light from an event 
is scattered by circumstellar or interstellar dust into the line 
of sight of the observer. LE imaging has long been used as a 
powerful tool for studying the three- dimensional dust struc- 
ture surrounding supernovae (SNe) dCrottsI 119881: IXu et al.l 
119951: ISugerman et al.ll2005t iKim et al.ll2008l) . In 1988, spec- 
tra of the inner and outer LE rings of SN 19 87A confirmed 
the re s emblance to a maximu m-light spectrum (Gouif fes et al.l 
119881: ISuntzeff et alJ fl988bh . More recently, however, this 
technique of targeted LE spectroscopy has been used on his- 
torical SNe to identify the spectral type of the original out- 
burst. After the serendipitous discovery of LEs from ancient 
SNe in the Large Mag ellanic Cloud (LM C) during the Su- 
perMACHO project by iRest et all (l2005al) . follow-up spec- 
troscopy by iRest et all (12008b!) identified the type of SN re- 
sponsible for the remnant SNR 0509-675. The spectral types 
of the Cas A and Tycho SNe have also be en identified using 
this method of targe ted LE spectro s copy (|Rest et ai1l2008al : 
iKrause et all l2008allbh . IRest et al.l (1201 lbl) used LE spec- 
troscopy to detect asymmetries in the outburst of the Cas 
A SN, while IRest et al.l (l2012bl) recently obtained LE spec- 
troscopy from the "Great Er uption" of r\ Carinae. We refer the 
reader to lRest et al.l d2012al) for a review of LE spectroscopy 
emphasizing these more recent results. Given the accelera- 
tion of the field of LE spectroscopy in recent years, a detailed 
study of spectra of the well-known LE system of SN 1987 A 
can provide a foundation for future studies, as well as provide 
new insight into the explosion of SN 1987 A. 

In the case of historical LEs, one has to use the photometric 
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and spectroscopic history of a differe nt SN to model the LE 
spectrum (see, e.g. IRest et al.ll201 lbl where the Cas A out- 
burst is modeled with the lightcurve and spectra of SN 1993 J 
and SN 2003bg). For SN 1987A, however, the exact spectral 
and photometric history of the LE source is known with high- 
precision, allowing observed LE spectra to be compared un- 
ambiguously to an isotropic scenario. In addition to acting as 
a test bed for LE spectroscopy theories, the near-circular LE 
rings of SN 1987 A allow the original outburst to be probed 
for asymmetries in an entirely direct way. Different position 
angles (PAs) on the LE system probe different viewing angles 
from which we can view the time-integrated spectrum of the 
original event. 

Observations as well as theoretical simulations indicate that 
core-collapse SNe are asymmetric in nature. Polarization 
measurements show deviations from sphe rical symmetry in 
all co re-collapse SN with sufficient data (IWang & Wheeler! 
120081) . while recent simulations in two and three spa- 
tial dimensions show large deviations from symmetry (e.g.. 
Hammer et ail 1201 Ot iGawrvszczak et al.l 120101: iMiiller et all 



20121) . Coupled with the fact that state of the art spheri- 



cally symmetric core-collapse simulations in one spatial di- 
mension fail_at_producing an explosion for > 10M Q pro- 
jectors (lBurrowsll2012l) . multi-dimensional p hysics such as 
the st anding accretion shock instability (SASL lBlondin et ail 
120031) may play a key role in understanding the explosion 
mechanism. As a new method for observing asymmetries, 
LE spectroscopy of SNe may be able to provide new insight 
into the origin of SN asymmetries and their relation to the 
core-collapse explosion mechanism. 

SN 1987 A was known to be an asymm etric Type IIP 
SN. Early p o larization measurem ents (e.g., iJeffervl 119871: 
iBailevI 119881: ICropper et al.l |1988|) and an elongated ini- 
tial speckle image (IPapaliolios et al.|[l989l) suggested a non- 
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spherical event. Fine structure in the Ha line (the "Bochum 
event") at 20-100 days after the explosion as well as red- 
shifted emission lines at more than 150 days after explo- 
sion provide strong evidence for radial-mixing of heavy el- 
ements into the upper envelope ([Hanuschik & Dachsl 119871: 
iPhillips & Heathcotel[T989t iSpvromilio et al.ll 19901) . Updated 
models of the bolometric lightcurve at early epochs also 
require radial- mixing of 56 Ni (IShigevama & Nomotol 1990): 
IUtrobinll2004l) . Direct HST imaging by IWang et al.l (12003 
showed an elongated remnant ejecta, claimed to be bipolar 
and showing evidence for a jet-induced ex plosion. Recent in- 
tegral field spectroscopy of the ejecta by iKiaer et all ( 120101) 
showed a prolate structure for the ejecta oriented in the plane 
of the equatorial ring, arguing against a jet-induced explosion. 

Here we present detailed imaging and spectroscopy of the 
SN 1987 A LE system. We infer the relative contributions of 
the different epochs of the SN to the observed L E spectrum 
by mod elling images of the LE, as demonstrated bv lRest et al] 
(1201 lal) . We fit the observed LE spectra and use specific ex- 
amples to illustrate the models ability to correctly interpret 
the observations. We then use the LE spectra to probe for 
asymmetries in the explosion of SN 1987 A and compare to 
previous observations of asymmetry. 

2. OBSERVATIONS AND REDUCTIONS 
2.1. Imaging 

Imaging of the SN 1987 A LE system was perfor med under 
the Su perMACHO Project microlensing survey dRest et alJ 
I2005bl) using the CTIO 4m Blanco telescope. The survey 
monitored the central portion of the LMC for five seasons be- 
ginning in 2001 using the 8K x 8K MOSAIC imager (plus 
atmospheric dispersion corrector) with the custom "VR" fil- 
ter (A c = 625 nm, <5A = 220 nm). Exposure times were be- 
tween 150 and 200 s. We have continued to monitor the 
field containing SN 1987 A and its LEs (field sm77) since 
the survey ended. Data reduction and difference images 
were perf ormed using the ESSENCE/SuperMACHO pipeline 
photp ipe dRest et al]|2005bt iGarg et alj[2007t iMiknaitis etafl 
120071) . A stacked and mosaiced difference image of the LE 
system is shown in Figure Q] The near-circular rings illumi- 
nate three general dust structures: a smaller structure ~ 85 pc 
in front of SN 1987 A, only visible in the south; a near- 
complete ring at ~ 185 pc and brightest in the north-east; 
as well as a larger and fainter near-circular illumination at 
<~ 400 pc in front of the SNR. All thre e of these struct ures 
have been mapped previously in detail bv lXu et ail (1 1995b and 
can be referred to as three "sheets" of ISM dust roughly in the 
plane of the sky. However, it is important to note that the LE 
flux observed in FigureQ]is due to dense filamentary structure 
within the general "sheets" of ISM dust. The physical proper- 
ties of the scattering dust filaments (inclination and thickness) 
can vary greatly within what appears to be a uniform "sheet" 
of dust. This distinction is important to make when modelling 
the LE photometry and spectroscopy. 

2.2. Gemini Spectroscopy 

We obtained multi-object optical spectroscopy (MOS) for 
five fields of the SN 1987 A echo system in the 2006B 
term, using the R400 grating and GG455 blocking filter on 
the Gemini Multi-Object Spectrograph (GMOS) on Gemini- 
South. SuperMACHO or GMOS preimages were differenced 
with previous SuperMACHO images to establish the precise 
location of the echo system, allowing the design of GMOS 




FIG. 1. — Mosaiced SuperMACHO difference image of August 2006 and 
January 2002 epochs of the SN 1987 A LE system. The locations of the 14 
GMOS slits appear in red, while the black cross indicates the location of SNR. 
Slitlets are not to scale. The black squares collect LEs which were observed 
in the same MOS field. North is up and east is to the left. 



masks. The locations of five MOS fields and the 14 1.0" wide 
LE slitlets are shown in Figure Q] 

The nod and shuffle mode of GMOS was used to best iso- 
late the diffuse LE signal from the nebular sky background. 
To ensure clear sky in the offset position, the telescope was 
nodded off-source several degrees every 225 s. On-source in- 
tegration times for the brightest echo fields were 30 minutes, 
and 90 minutes for the fainter, more diffuse echoes. Dur- 
ing each night of observations, flat field and CuAr spectral 
calibration images were taken before or after the science im- 
ages. We also acquired dark images with the same nod and 
snuffle parameters used during integration. These specialized 
dark images allow features associated with charge-traps in the 
GMOS CCDs to be masked out during reduction. 

The observations result in LE spectra from SN 1987 A with 
a spectral range of 4500-8500 A, contaminated by strong 
emission lines from the nebular region, as well as any stellar 
continuum entering the slitlets. The P-Cygni Ha line seen in 
the spectra is further degraded due to six strong nebular emis- 
sion lines ([Nil] A6548, Ha A6563, [Nil] A6585, Hel A6678, 
[Sill] A6717, [Sill] A6731). Signal from the emission peak 
of the P-Cygni Ha line is therefore lost to nebular emission, 
and this region of the profile cannot be interpreted with con- 
fidence. In particular, the location and flux of the emission 
peak cannot be determined without interpolations which we 
found introduced large uncertainties. 

The nebular contamination can be greatly reduced since the 
LE system is expanding superluminally on the sky (as seen 
in Figure [TJ. We therefore have the unique opportunity in 
astronomy to observe only the background sky signal at the 
same location on the sky as the initial sky+object observa- 
tion. Sky-only observations were taken in the 2009A and 
2009B semesters for the five MOS fields. The spectroscopy 
was carried out using the same GMOS configuration as the 
previous sky+object 2006B observations, with 30 minutes on- 
source integration times in nod and shuffle mode. The spec- 
trophotometric standard star LTT 3864 was observed during 
the 2009A observations. 
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2.3. Spectra Reduction 

All spectra were reduced using IRAF0 and the Gemini 
IRAF package, with bias subtraction and trimming done on 
each CCD individually. Science images were dark subtracted 
using the special nod and shuffle darks, with each CCD be- 
ing cleaned of cosmic rays using the Lapl acian rejection al- 
gorithm LACOSMIC (Ivan Dokkuml 1200 ll) . Nod and shuf- 
fle sky subtraction was performed on each science exposure 
using the gnsskysub task. Individual slitlets were automati- 
cally cut from the wavelength calibrations frames, however 
this procedure required manual tweaking to achieve the best 
cutting possible. A wavelength solution was found for each 
slitlet in the calibration frame. Science frames were flattened, 
CCDs mosaiced together, and slitlets cut from the images. 
The wavelength solution for each slitlet was then applied to 
the corresponding science spectrum in two dimensions. Since 
the LE signal is so weak, determining an extraction aperture is 
very difficult. Instead, we collapsed the entirety of each slitlet 
to one dimension using a block average. 

No suitable spectrophotometric standard star was observed 
during the 2006B observing semester, nor did a suitable stan- 
dard star matching our unique instrument configuration exist 
in the data archive nearby temporally. However, the spec- 
trophotometric standard LTT 3864 was observed during the 
2009A observations. We therefore used the 2009A standard 
to perform the flux calibration for all of the 2006 and 2009 
spectra. While this leads to large errors in absolute flux cal- 
ibration, our analysis is independent of absolute flux levels. 
Instead, the purpose of the flux calibration used here is to re- 
move the effect of CCD sensitivity from the spectra, which is 
relatively stable in time. A LE spectrum is a weighted aver- 
age of many epochs from the outburst. Determining the con- 
tinuum of such a spectrum is therefore difficult, and the alter- 
native analysis procedure of removing the continuum from all 
of our spectra would have had a larger error than that of the 
flux calibration. 

The sky-only spectra from 2009 were subtracted in one di- 
mension from the 2006 LE spectra for each MOS slit. This 
procedure was performed interactively, adjusting the scale and 
wavelength of the sky spectrum in small increments using 
the IRAF task skytweak to minimize sky subtraction residu- 
als in the final spectra. Since we focus our science analysis on 
the Ha line, we performed the sky subtraction such that the 
subtraction residuals in the P-Cygni profile of Ha were mini- 
mized. All spectra were Doppler-corrected, adopting a LMC 
radial velocity of 286.5 km s" 1 dMeaburn et alJI 19951) . 

3. ANALYSIS 

3.1. Light Echo Spectra 

The upper panel of Figure [2] shows a reduced LE spectrum 
from the 2006 observations along with the underlying sky- 
only nebular spectrum from 2009. The high-velocity P-Cygni 
profiles including Ha resemble a Type II SN spectrum at max- 
imum light, confirming the LE spectra probe the original out- 
burst light of SN 1987 A. Strong, narrow nebular lines not as- 
sociated with the SN outburst dominate the spectrum around 
5000A and the Ha line at 6563 A. The result of performing 
difference spectroscopy is shown in the lower panel of Fig- 
ure|2] where the sky-only signal has been subtracted from the 

4 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 
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FIG. 2. — Upper Panel: Example LE+sky spectrum from 2006 shown 
with the corresponding sky-only nebular spectrum from 2009. Lower Panel: 
The result of difference spectroscopy, showing the LE-only spectrum along 
with remaining nebular residuals from the subtraction. The strength of the 
Hce emission peak can now be estimated in the difference spectrum. 

2006 echo spectrum. Flux from the emission component of 
the Ha line is recovered that was initially lost to nebular con- 
tamination. In the majority of our spectra we are able to re- 
solve the Ha emission peak without loss due to the nebular 
lines. 

3.2. Modeled Isotropic Spectrum 

To search for asymmetries in SN 1987 A using spectroscopy 
of its LEs, we need to first consid er what def i nes asy mmetry 
in observed spectra. As shown in iRest et aD (1201 lal) . multi- 
ple LE spectra from the same transient source, with differing 
line strength ratios, does not necessarily imply an asymmetry 
in the outburst light. An observed LE spectrum will depend 
strongly on the physical properties of the scattering dust as 
well as the instrument configuration and seeing conditions at 
time of observation. Specifically, the dominant epochs of the 
original outburst probed by a LE spectrum can change with 
the above properties. Relative comparisons of LE spectra can 
therefore easily lead to false-detections of asymmetry if the 
above factors are not taken into consideration. Each LE spec- 
trum should instead be compared to a single isotropic spec- 
trum modelled for each observation. We define our isotropic 
model spectrum as the original outburst of SN 1987 A, as it 
would have been spectroscopically observed if scattered by 
dust corresponding to the observed LE. This provides a di- 
rect probe of asymmetry since differences between observed 
and model spectra represent deviations from the historically 
observed outburst of SN 1987 A. 

A LE observed on the sky has a flux profile (flux versus 
p, where p is the distance on the sky from the SNR, see up- 
per right panel of Figure |3]l that is the projected lightcurve 
of the source event stretched or compressed depending on the 
inclination of the scattering dust. The width of the scattering 
dust and the point spread function (PSF) further distort the 
observed LE. The dust inclination is measured through imag- 
ing the LE at multiple epochs and can be measured prior to 
spectroscopy. The PSF is also known, allowing the width of 
the scattering dust filament to be determined through fitting 
the LE flux profile. 

The inclination, a, is obtained by monitoring the apparent 



4 



Sinnott et al. 



Projected Phase in days 
,800 400 -400 -800 





472 473 474 475 476 

z in ly 



165 170 175 

p in arcsec 



1.0 
0.8 

"|> 0.6 

ii 
S 

U 0.4 

0.2 
0.0 








- - Original LC 
Effective LC 





-100 100 200 300 400 
Phase in days 



200 400 600 
Phase in days 




FIG. 3. — Summary of modeling procedure for each LE spectrum. The 
inclination of the dust filament is measured (upper left panel), allowing the 
shape of the LE profile to be modeled to determine the dust width (upper 
right panel). The properties of the slit determine the window function con- 
taining the relative contribution from each epoch (middle left panel), which 
when multiplied by the original lightcurve gives the effective lightcurve for 
this observation (middle right panel). We integrate the historical spectra of 
SN 1987A, flux- weighted with the effective lightcurve, and take into account 
wavelength-dependent scattering and reddening to obtain a fit between the LE 
model and the observed spectrum (lower panel). The lower panel plot shows 
the model LE fit (solid red) along with the attempted fit with the maximum 
light spectrum of SN 1987A (dashed blue) and the full lightcurve-integrated 
spectrum (dotted black). 

motion of the LE on the sky. Positive values of a correspond 
to scattering dust sheets tilted out of the plane of the sky, 
from the positive p axis toward the negative z axis, where z 
is the distance in front of the remnant. Using the SuperMA- 
CHO database of difference images around SN 1987 A, we are 
able to monitor the apparent motion of the echo system over 
seven years from 2002 to 2008. For each difference image 
and each slit location in FigureQ] we measure the one dimen- 
sional profile of the difference flux as a function of distance 
away from the SNR, p. Carefully monitoring the superlumi- 
nal apparent motion in this way allows us to determine the 
inclination of the scat t ering d ust for each slit location, as de- 
scribed in iRest et alj d201 lab (see, e.g., upper left panel of 
Figure |3). We find the apparent motion (and therefore incli- 
nation) is typically stable over a one year period. We find 
typical uncertaint ies in the inclinatio n are less than 5° for our 
14 dust locations. IRest et ail d201 lal) showed that, for the case 
of SN 1987 A, changes in inclination on the order of 10° do 
not alter the final modeled spectrum in a significant way. With 
the inclination known, the photometric history of SN 1987 A 
dHamuv et alJll988HSuntzeff et al.lll988atlHamuv & Suntzeffj 



fl990h can be used to model the one dimensional flux profile 
and determine the best-fitting scattering dust width. The up- 
per right panel of Figure [3] shows an example fit to the LE 
profile. 

The above procedure determines the properties of the scat- 
tering dust for each LE location. We also take into account 
the location, orientation, and size of the spectroscopic slit to 
compute a window function that is unique to each LE obser- 
vation (see, e.g., middle left panel of Figure [3]). This win- 
dow function is the relative contribution from each epoch of 
the lightcurve when observed through the slit. The slit offset, 
A/5, measures the offset between the peak of the LE profile 
and the location of the slit on the sky. This is an important 
parameter, as it shifts the window function in time, systemati- 
cally probing later or early epochs in the spectrum if the slit is 
positioned closer or further away from the SNR, respectively. 

Multiplying the window function with the original 
lightcurve of SN 1987 A, we determine an effective lightcurve 
corresponding to each observed LE spectrum (see, e.g., mid- 
dle right panel of Figure [3]). The isotropic model spec- 
trum is then the integration of the historical spectra of 
SN 1987A weig hted with t he effe ctive lightcurve. We re- 
fer the reader to IRest et aD (1201 lal) where our model is de- 
scribed in detail with SN 1987A examples. To integrate 
the historical spectra of SN 1987A, we use the integration 
method described in IRest et all d2008bl) . wit h the spectral 
datab a se from SAAO and CT I O obs e rvations dMenzies et alj 
19871: iCatchpole et all fl987l fl9M iWhitelock et alJ 119881: 



Catchpole et al.lll989HWhitelock et alJ l!989t ICaldwell et al] 
1993tlPhillips et alJI 198811 1990h . 

The wavelength-dependent effect of scattering by dust 
grains is taken into account when fitting the isotropic model 
to our observed LE spectra. We calculate the integrated scat- 
tering function for each LE (i . e. scat tering angle) using the 
method described in lRest et"afl d2012al), using the "LMC avg" 
carbon aceous-silicate grain model of IWeingartner & Drainel 
( 120011) . We then fit the model spectrum to the observed LE 
spectrum (see, e.g., lower panel of Figure[3]) varying three pa- 
rameters: the reddening (E(B-V) assuming R v = 3.2), a nor- 
malization constant, and a small flux offset (to account for 
small errors in sky subtraction). 

Figure [3] summarizes the m odeling procedure outlined 
above and in lRest et al.l d201 lal) that occurs for each LE spec- 
trum. The lower panel shows the final isotropic model spec- 
trum in red, which fits the features in the observed LE spec- 
trum very well. To highlight the importance of the model- 
ing procedure, we also attempt to fit the LE spectrum with 
a spectrum of SN 1987 A near maximum light (dashed blue 
line), and an integrated spectrum obtained by flux-weighting 
with the original lightcurve of SN 1987A (dotted black line). 
The maximum light spectrum cannot reproduce the strength 
of the observed Ha emission. Conversely, the full lightcurve- 
integrated spectrum has an excess in Ha emission and a lower 
Ha velocity in the absorption minimum. Only by integrating 
using the modeled window function (solid red line) are we 
able to obtain a good fit to the Ha profile. Below we high- 
light the importance of the LE modeling by considering two 
scenarios. 

3.2.1. Dust-Dominated Scenario 

Both the physical properties of the scattering dust as well 
as the configuration of the observation affect the observed LE 
spectrum. Here we highlight the effect of the scattering dust. 
Table Q] shows the observed inclinations and best-fitting dust 
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widths for the 14 LE locations shown in Figure Q] Figure [4] 
shows two examples of observed flux profiles, corresponding 
to the two LEs in FigureQ]at PA ~ 1 15°, LEI 13 and LEI 17. 
Note the names of the LEs correspond to the PA of the LE 
with respect to the SNR. Both LEs are at essentially the same 
PA along the echo system, so probe the same viewing angle 
onto the SN. However, as shown in Figure[T]the two slits are 
placed on physically distinct dust filaments. The bright fila- 
ment in LEI 13 is 53 ly closer to the observer along the line 
of sight than the bright LEI 17 filament, while both are 3 ly 
thick. The observed widths and shapes of the LE profiles in 
Figured] are also different, with LEI 13 having a larger width 
on the sky and a more symmetric shape. This is explained 
entirely by the ~ 50° difference in the inclinations of the two 
dust filaments. The highly inclined dust filament for LEI 13 is 
within ~ 15° of the tangential to the scattering ellipsoid and 
compresses the lightcurve on the sky, causing the LEI 13 pro- 
file to resemble a point-like source with a broad, symmetric 
profile. The lower inclination of the LEI 17 dust filament pre- 
serves the shape of the lightcurve with a broad increase in flux 
at small p values (the long decay of the SN lightcurve) and a 
sharp decrease in flux at large p values (the short rise of the 
SN lightcurve). It should be noted that the LEI 13 profile in 
the left panel of Figure H shows dust substructure to the left 
of the main LE peak. The additional contribution from the 
secondary peak is weak pre-maximum flux and we found the 
effect on the Ha emission strength in the integrated spectra 
to be < 4%. However, we do take substructure in the dust 
into accou nt in the extreme case of LE186 as shown in Sec- 
tion l3.3.l1 below. The effect of dust substructure is considered 
more formally in the Appendix. 

The modeled isotropic spectra are shown in Figure [5J with 
the upper left panel showing the observed LE spectrum for 
the LEI 17 slit (black), the corresponding LEI 17 model (thin 
red), and the non-matching LEI 13 model (thick cyan) de- 
rived using the larger inclination. Since the dust inclination 
of LEI 13 compresses the lightcurve on the sky, a much larger 
range of epochs span the width of the slit resulting in a wide 
window function compared to the LEI 17 slit (Table Q] lists 
the approximate range of epochs probed by each LE). The 
LEI 13 model therefore has an excess of Ha emission from 
the inclusion of late-time nebular epochs. The LEI 17 model 
is able to correctly match the observed LEI 17 Ha strength. 
The lower panel of Figure [5] shows the same LEI 13 model 
with the corresponding LEI 13 observed LE spectrum, show- 
ing good agreement. The fact that two distinct LEs with dif- 
fering line strengths and differing dust properties show the 
same result - that the observed LE spectra from this viewing 
angle can be modeled with an isotropic historical spectrum 
- shows the int erpretation of LE s pectroscopy described here 
and in detail in lRest et alj ( 1201 lal) is correct. 

3.2.2. Obseri'ation-Dominated Scenario 

Here we consider a scenario where the properties of the ob- 
servation (specifically the slit location) are the dominant fac- 
tors in two very different observed LE spectra that have sim- 
ilar viewing angles onto the photosphere. Figure [6] shows the 
LE profiles of LE053 and LE066, corresponding to the slits 
at PA^ 60° in FigureQ] The gray shaded region indicates the 
location of the slit for both profiles. The slit of LE066 was 
placed ~ 1" farther away from the SNR than the peak of the 
LE. This was unintentional since the location of the LE peak 
may not be very well known when making the MOS masks. 
The offset results in very different observed LE spectra for the 
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FIG. 4. — Observed LE profiles and corresponding best-fit models (red line) 
for LEI 13 (left) and LEI 17 (right). The larger width and symmetric shape of 
the LEI 13 profile on the sky is due to the dust filament being inclined ~ 50° 
closer to the tangential to the line-of-sight than the LEI 17 dust filament, com- 
pressing the entire lightcurve on the sky and resulting in a more symmetric 
point-like LE profile. 
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FIG. 5. — The upper left panel shows the observed LE spectrum of LEI 17 
with its corresponding isotropic model, along with the incorrect model of the 
nearby dust filament LEI 13. Only the matching LEI 17 model can fit the 
observed LE spectrum, while the lower left panel shows the LEI 13 observed 
LE spectrum can be fit with its corresponding model. The right panels show 
enlargements of the Ha profile in both cases. 

two profiles, as shown in Figure|7] When compared to the ob- 
served LE spectra of Figure|5]and LE053 in Figure|7] the spec- 
trum of LE066 looks more like an early-time spectrum with a 
lower Ha emission-to-absorption ratio and a higher tempera- 
ture continuum. Since the slit is centered on the pre-maximum 
light portion of the LE profile, LE066 does not probe epochs 
> 100 days after maximum light, while LE053 probes epochs 
out to ~ 300 days after maximum light. As in the previous 
example, both model isotropic spectra are good fits to their 
observed LE spectra. However, the model of LE053 has a 
lower Ha absorption velocity and a larger Ha emission-to- 
absorption ratio, making it a poor fit to the LE066 observed 
LE spectrum. 

The above dust- and observation-dominated examples high- 
light the danger in comparing LE spectra without suffi- 
cient modeling. Without the corresponding models for the 
LE053 and LE066 observed spectra, one might make a false- 
detection of asymmetry. Specifically, that the LE066 spec- 
trum has a larger Ha velocity and excess in emission com- 
pared to the LE053 spectrum. If the two LEs occurred at op- 
posite PAs this asymmetry claim would be even more tempt- 
ing to make. 

3.3. Evidence for Asymmetry in SN 1987A 

Each PA on the LE system of SN 1987 A represents a dis- 
tinct viewing angle with which to view the original outburst. 
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TABLE 1 



LEID a 


PA b 


Zo c 


e i 


a c 




Ap ffset B 


Epochs' 1 




(°) 


(ly) 


(°) 


(°) 


(ly) 


(") 


(days) 


LEO 1 6 


15.8 


428.91 


17.1 


22 ± 6 


4.4y ± U.UJ 


-0.35 


n 1 on 
0- 1 oU 


LE029 


29.2 


473.13 


16.3 


28 ± 2 


a a c\ i n m 

4.49 ± 0.03 


-0.19 


5-200 


LE032 


31.8 


483.02 


16.2 


-46 ± 2 


o T 1 I rt etc 

2.71 ± 0.05 


-0.27 


30-180 


LE034 


33.9 


486.53 


It). 1 


-lo ± z 


i in _L n i ^ 
3.3 / ± 0. 1 J 


+0.26 


0-140 


LE053 


53.0 


505.05 


15.8 


38 ±7 


5.26 ± 0.03 


-0.18 


0-215 


LE066 


66.3 


590.39 


14.7 


11 ±2 


4.50 ± 0.55 


+1.05 


0-80 


LE069 


69.3 


624.64 


14.3 


66.8 ± 0.2 


1.80 ± 0.07 


+0.57 


0-110 


LE076 


76.4 


648.50 


14.0 


11 ±4 


0.60 ± 0.22 


-0.08 


40-135 


LEI 13 


112.5 


671.04 


13.8 


79 ±1 


3.06 ± 0.09 


-0.39 


0-330 


LEI 17 


116.8 


618.40 


14.3 


33 ±5 


3.10 ±0.39 


-0.14 


35-155 


LEI 80 


180.2 


289.69 


20.7 


-49 ±2 


4.53 ± 0.06 


+0.73 


0-290 


LEI 86 


185.6 


270.34 


21.3 


70 ±3 


3.52 ±0.18 


-1.00 


0-420 


LE325 


325.3 


369.67 


18.4 


-6± 11 


1.13 ± 0.16 


+0.18 


20-120 


LE326 


326.1 


370.33 


18.4 


15 ±3 


2.24 ±0.51 


+0.21 


0-130 



d Naming convention corresponds to position angle of LE slit with respect to SNR. 

Position angle of the scattering dust with respect to the SNR. 
G Distance along line of sight from SNR to scattering dust. 
^ Scattering angle. 

c Inclination of scattering dust with respect to the plane of the sky. a increases from the positive p axis 
towards the negative z axis. 



Average best-fitting width of scattering dust sheet over the length of the spectroscopic slit. 



g Offset between center of slit and peak LE flux, where positive values indicate slit further from SNR than 
peak echo flux. 

k Approximate range of epochs probed by the LE. Corresponds to epochs in window function in which the 
relative contribution contribution is > 50% (prior to flux-weighting) in the LE integration. 




FIG. 6. — Observed LE profiles and corresponding best-fit models (red line) 
for LE053 (left) and LE066 (right). The grey shaded regions indicate the 
location of the spectroscopic slit. The slit offset from the LE peak in the 
LE066 LE produces an observed spectrum that is heavily weighted towards 
the earlier epochs of the outburst. 

LEs originating to the north of SN 1987 A probe outburst 
light originating from the northern portion of the photosphere. 
Therefore, by observing the LEs as a function of increas- 
ing PA, we are able to view the outburst of SN 1987 A as a 
function of north to south lines of sight. The opening angles 
probed by the LE ring are 28° -43°, with the scattering angles 
listed in Table Q] 

Observed LE spectra and corresponding dust-modelled 
isotropic spectra are plotted as a function of PA in Figures [8} 
fl2l with each figure representing a field as shown in Figure Q] 
Figure [8] corresponds to the most northern field. As previ- 
ously mentioned, we should avoid searching for differences 
between observed LE spectra. Instead, deviations in asymme- 
try are defined by deviations from the dust-modeled isotropic 
spectrum (red line) for each line of sight. A deviation from 
the model then describes a deviation from the appropriate 
weighted set of outburst spectra of SN 1987 A as it was histor- 
ically observed along the direct line of sight. 

3.3.1. Ha Profiles 

Figures[8lfl2lshow that most features and line strengths ob- 
served in the optical LE spectra can be fit with the isotropic 
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FIG. 7. — The upper left panel shows the observed LE spectrum of LE066 
with its corresponding isotropic model, along with the incorrect model of the 
nearby dust filament LE053. Only the higher velocity LE066 early-epoch 
model is a good fit for the observed LE spectrum. The lower left panel shows 
the LE053 observed LE spectrum which probes later epochs of the explo- 
sion. The lower velocity and larger Ha emission-to-absorption ratio of the 
LE053 model is a good fit to the observed spectrum. The right panels show 
enlargements of the Ha profile in both cases. 



historical spectrum of SN 1987A, without the need to invoke 
asymmetry in the outburst. However, the fine-structure and 
strength of the Ha A6563 line shows deviations from symme- 
try as a function of PA (i.e. viewing angle onto the SN). 

Figure Qj] shows a closeup of the Ha profiles of seven 
unique viewing angles as a function of increasing PA from 
top to bottom. The most northern line of sight LE, LEO 16, 
shows a strong blue knee in the profile at ~ -2000 km s" 1 in 
addition to a strong excess of emission that is redshifted 
from the rest wavelength by ~ +800 km s" 1 (and by <~ 
+500 km s" 1 compared to the historical spectrum). The most 
southern LE, LEI 86, has a line of sight almost directly op- 
posite that of LE016. Its Ha profile (second from bottom 
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FIG. 8. — Observed LE spectra with corresponding modelled isotropic spec- 
tra for field at PA ~ 30° in Figure[T] Spectra are plotted from top to bottom as 
a function of increasing position angle as shown in Figure[T] LE016, LE029, 
LE032, LE034. Gaps in the observed spectra correspond to areas contami- 
nated by sky-subtraction residuals. 
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FIG. 9. — Observed LE spectra with corresponding modelled isotropic spec- 
tra for field at PA ~ 65° in Figure[T] Spectra are plotted from top to bottom as 
a function of increasing position angle as shown in Figure[T] LE053, LE066, 
LE069, LE076. Gaps in the observed spectra correspond to areas contami- 
nated by sky-subtraction residuals. 

in Figure [PJt shows a red knee in the fine-structure at ~ 
+2200 km s" 1 and an excess in emission shifted towards the 
blue by ~ -500 km s -1 (and by ~ -1000 km s" 1 compared to 
the historical spectrum). The qualitative difference in the 
Ha profile shape between the two extreme viewing angles, 
LEO 16 and LEI 86, is shown in Figure [14] This figure does 




FIG. 10. — Observed LE spectra with corresponding modelled isotropic 
spectra for field at PA ~ 115° in Figure [T] Spectra are plotted from top 
to bottom as a function of increasing position angle as shown in Figure [T] 
LEI 13, LEI 17. Gaps in the observed spectra correspond to areas contami- 
nated by sky-subtraction residuals. 
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FIG. 11. — Observed LE spectra with corresponding modelled isotropic 
spectra for field at PA ~ 185° in Figure [T] Spectra are plotted from top 
to bottom as a function of increasing position angle as shown in Figure [T] 
LE180, LE186. Gaps in the observed spectra correspond to areas contam- 
inated by sky-subtraction residuals. LE180 (upper spectrum) represents an 
example where our LE fitting algorithms fail to extract meaningful results 
due to the complexity and low signal-to-noise of both the spectrum and the 
LE flux profile. 



Observed LE Spectrum 
Modelled Isotropic Spectrum 




FIG. 12. — Observed LE spectra with corresponding modelled isotropic 
spectra for field at PA ~ 325° in Figure [T] Spectra are plotted from top 
to bottom as a function of increasing position angle as shown in Figure [T] 
LE325, LE326. Gaps in the observed spectra correspond to areas contami- 
nated by sky-subtraction residuals. 
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FIG. 13. — Closeup of observed Ha profiles (black) and corresponding 
dust-modelled isotropic spectra (red), as a function of PA for seven view- 
ing angles. The top Ha profile corresponds to the most northern slit in 
the upper right legend, and the profiles increase in PA from top to bottom: 
LE016, LE029, LE053, LE076, LEI 13, LE186, LE326. Residuals between 
the observed LE spectra and the dust-modeled isotropic spectra are plotted 
in blue below each profile. To prevent guiding the eye, the portions of the 
Ha profile containing residuals from sky subtraction have been masked out. 
The most northern profile, LEO 16, has an excess in redshifted Ha emission 
and a blue knee at ~ —2000 km s~' are observed. At almost opposite view- 
ing angle (LEO 186, second from bottom), the opposite asymmetry is ob- 
served in the Ha profile: an excess in blueshifted emission and a red knee at 
~ +2200 km s _I . Both asymmetry features diminish as afunction of viewing 
angle away from the north. Portions of the observed Ha peak contaminated 
by sky-subtraction residuals have been removed to prevent guiding the eye. 
No smoothing has been applied to any of the spectra. 




Radial Velocity {x in 3 km s '] 



FIG. 14. — Observed Ha lines from LE016 and LE186. Emission peaks 
have been interpolated with high-order polynomials. Spectra are scaled and 
offset for comparison purposes, as well as smoothed with a boxcar of 3 pixels. 
Although this plot does not take into account the important differences in LE 
time-integrations between the spectra, it highlights the overall difference in 
fine-structure in the two LE spectra from opposite viewing angles. Observing 
Ha profiles with opposite asymmetry structure at opposite viewing angles is 
surprising considering the opening angle between the two LEs is < 40° . 




Radial Velocity ( xlrj :i km s- 1 ) 



FIG. 15.— Historical CTIO Ha profile from IPhillips et alj f[988l) of 
SN 1987A taken on March 29, 1987, 34 days after explosion. The fine- 
structure, showing blue and red features at ~ ±3500 km s , defines the 
"Bochum event" which was observed 20-100 days after explosion. Evidence 
for this fine structure is not present after the histo rica l spectra are integrated 
with the LE scattering model (red l ines in Figure [\3\ . Ha fine structure ob- 
served in the LE spectra in Figure [T3l must therefore require stronger line 
profile asymmetry than the original "Bochum event" or the asymmetry must 
last for much longer than originally observed for SN 1987A. 

not take into account the effects of the LE observations on 
the integrated spectra and so cannot be used as a direct com- 
parison between the viewing angles. However, it does quali- 
tatively describe the asymmetry in the fine-structure between 
the two viewing angles. 

The two asymmetries in the Ha profile (summarized in Fig- 
ure n~4T> are the fine-structure, with a red knee in the north- 
ern profiles and a blue knee in the southern profile, as well 
as an excess of Ha emission that is redshifted in the north 
and blueshifted in the south. Both of these asymmetries ap- 
pear to be physical for a number of reasons. The excess in 
Ha emission as well as the blue knee feature smoothly dimin- 
ish in Figure[T3]as PA is increased. At PA <~ 1 10° in the south- 
east quadrant, both the Ha strength and fine-structure from 
LEs LEI 13 and LEI 17 can be fit with the historical model 
spectrum of SN 1987A. The fact that both Ha asymmetries 
(fine-structure and redshifted emission) smoothly diminish as 
a function of viewing angle from north to south in the eastern 
half of the LE ring is strong evidence that these asymmetries 
are physical. 

This fine-structure in the Ha profile is similar to the 
"Bochum event" o r iginal ly observed in SN 1987A 
dHanuschik & Dachsl fl987t IPhillips & Heathcotd [19891) . 
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FIG. 16. — LEFT: Observed LE profile on the sky for the southern LE186. 
Over-plotted in red are models for two scattering dust filaments which both 
contribute signal to the slit (plotted as the shaded region). Although the nar- 
row filament at larger p contributes late-time epochs to the spectrum, al- 
lowing for a larger emissio n-to -absorption ratio, the excess emission in the 
LEI 86 Ha profile in Figure [T3l cannot be modeled with the historical spectra 
of SN 1987A. See Appendix for full consideration. RIGHT: Observed LE 
profile of LEI 80. See Appendix for full consideration on why the LE fitting 
algorithm does not provide meaningful results in this complex case. 

where blue and red satellite emission features were ob- 
served 20-100 days after explosion. Figure [151 shows the 
Ha profile of a CTIO spectrum taken 34 days after explosion 
dPhillips et aljfT988l) . The blue and red emission satellites are 
similar to the blue and red knees observed in the LEO 16 and 
LEI 86 LE spectra, respectively. It is tempting to compare 
the Ha fine-structure observed in Figure Q3] directly with 
the fine-structure observed in the LE profiles. However, LE 
spectra represent an integration of many epochs of signal 
and cannot be compared directly to a single epoch spectrum. 
Although the "Bochum event" is prominent in Figure Q3] 
there is no evidence for fine-structure in any of the isotropic 
model spectra in Figure [13] That is, the "Bochum event" 
does not survive the smoothing effect when the historical 
spectra of SN 1987 A are integrated. This can be seen in 
the lower panel ofFigure [3] where the fine-structure from 
the "Bochum event" is faintly visible in the maximum 
light spectrum, but entirely nonexistent in the integrated 
model spectra. The fact that fine-structure similar to the 
"Bochum event" is present in the observed LE spectra is 
evidence for underlying fine-structure that is much stronger 
than that originally observed for SN 1987 A. In general, 
small deviations from symmetry observed in a LE spectrum 
must represent very large deviations from symmetry in the 
underlying outburst in order to remain after the integrating 
effect of the LE phenomenon. This is an important aspect of 
all LE spectroscopy work that is rarely emphasized. 

The excess emission in the north and south from profiles 
LE016 and LE186, respectively, cannot be explained within 
uncertainties. To obtain the largest temporal coverage we used 
both CTIO and SAAO spectra when integrating the model 
spectrum. There exists known systematic differences between 
the d atasets of the origin al observations from the two loca- 
tions (lHamuy et al.l ll990b. Using the two data sets individ- 
ually in our analysis pipeline can lead to differences in the 
Ha strength by 5%- 10% in our models, but cannot account 
for the 30% -40% excess observed in the LE016 and LE186 
spectra. We also stress that our analysis is based on relative 
flux comparisons only. 

The LE profile on the sky for LEO 16 is very similar to that 
shown in Figure [3] Since it is a single peak, the theoreti- 
cal maximum amount of Ha emission in the model would 
correspond to integrating the historical spectra with the full 
lightcurve of SN 1987A out to t = t now , rather than an effective 
lightcurve that is truncated by a window function (i.e. having 
an infinitely thick dust sheet). However, even this unphysical 
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FIG. 17. — Observed Fe II A5018 line in LE016 most-northern LE spec- 
trum. The isotropic dust-modelled spectrum is over-plotted (solid red line), 
along with the same isotropic spectrum shifted 800 km s to the red (dashed 
blue line), corresponding to the best fit to the peak location of the A5018 
line and of similar magnitude to the velocity shift observed in the Ha profile. 
Note that the velocity-shifted model spectrum is a better fit to the observed 
LE spectrum, indicating that the Ha velocity asymmetry may be correlated 
with the velocity of the Fe II A5018 line. 

limit cannot reproduce the Ha emission-to-absorption ratio 
that is seen in the observed LE profile of LE016. The LE pro- 
file on the sky and the slit location for LEI 86, which has two 
closely spaced filaments, are shown in Figure [T6] Since the 
slit was placed on the first peak, there will be late-time neb- 
ular emission entering the slit from the second peak at larger 
p. In such a case, it is possible to obtain a larger emission- 
to-absorption ratio. However, we have taken this into account 
by modeling both LE peaks and determining the relative con- 
tributions from each peak entering the slits. This effect, dis- 
cussed in more detail in the Appendix, cannot account for the 
excess in emission in LEI 86. We also stress that the apparent 
lack of any Ha emission excess in LEI 80 is due to limitations 
in our LE fitting algorithrims when dealing with such low 
signal-to-noise data. LEI 80 is discussed in more detail in the 
Appendix. It should be noted, however, that LEI 80 appears 
to show the same fine-structure in Ha (blueshifted emission 
peak and red knee) despite the low signal of the spectrum. 

The opening angles probed by the LE lines of sight are < 
45° . The fact that large asymmetries are seen in the observed 
integrated LE spectra for relatively small changes in viewing 
angle is evidence for a strong asymmetry in the explosion. 

3.3.2. Additional Evidence of Asymmetry 

Although the Ha profile shows the strongest asymmetry 
signature that is dependent on viewing angle, we note here 
two additional possible sources of asymmetry from the LE 
spectra: a velocity shift in the Fe II A5018 line that appears 
to be correlated with the Ha asymmetry, and an unidentified 
feature near 5265 A present in only one direction. 

Figure [TT] plots the Fe II A5018 line for the observed LE 
spectra of LEO 16 (black), the asymmetric northern viewing 
angle. The dust-modelled isotropic spectrum is plotted as 
solid red. We also plot in dashed blue the same isotropic spec- 
trum redshifted by +800 km s -1 , which gives the best fit to the 
emission peak of the Fe II A5018 line (determined by eye). 
The redshifted spectrum is a better fit to the line and is consis- 
tent with the 600-800 km s" 1 redshift observed in the Ha line 
for the same viewing angle. As with the Ha line asymmetry, 
the magnitude of the best- fitting redshift for the Fe II A5018 
line decreases as a function of viewing angle away from the 
LEO 16 line of sight. For the case of LEI 86, the line is weak 
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FIG. 18. — Observed LE spectrum for LE016 (top) compared with the 
two "equatorial" LE locations (LEI 13 and LEI 17, middle and lower, respec- 
tively) shown at PA ~ 115° in Figure [T] Both spectra show an unidentified 
feature near 5265A that is not present in any other viewing angle. The fact 
that the line appears broader in the upper spectrum is consistent with that 
LE probing a narrower range of early, high-v elocity epochs compared to the 
lower spectrum as described in Section l3.2.1l 

and heavily contaminated from sky residuals. Although it has 
the most blueward Fe II A5018 line of all observed LE spec- 
tra, it is unclear if the line is actually blueshifted with respect 
to the isotropic model or to zero. The structure is complicated 
since this line is part of a blend of the Fe II AA4924, 5018, 
5169 features. 

Figure fT~8l shows an unidentified feature near 5265 A in the 
observed LE spectra for both "equatorial" (i.e. perpendicular 
to maximum asymmetry axes 16° /186° ) lines of sight at PA 
~ 1 15° . Many LE spectra have signal-to-noise ratios higher 
than the lower spectrum of Figure[T8] however there is no ev- 
idence for the feature in any other line of sight. The fact that 
the feature is broader in the LEI 17 (lower) spectrum is con- 
sistent with that LE probing a narrower range o f early, high- 
velocity epochs as described in Section 13.2.11 Considering 
this wavelength region is populated with many line blends and 
the fact that the spectra represent a time-integration, it is cur- 
rently unclear if the feature is due to differences in chemical 
or velocity properties at this viewing angle. 

4. DISCUSSION 

Ignoring the smaller knee-like fine-structure, the 
Ha profiles of the LEO 16 and LEI 86 viewing angles ar- 
gue strongly for a one-sided asymmetry in SN 1987 A. An 
overabundance of 56 Ni in the southern far hemisphere would 
create an excess in nonthermal excitation of hydrogen. 
This results in an excess in redshifted emission for the 
northern LE016 viewing angle. If the overabundance of 
56 Ni is inclined close to the plane of the sky (within 21°), 
the overexcitation would occur in the near hemisphere with 
respect to the LEI 86 line of sight, explaining the blueshifted 
Ha emission in LEI 86. The question remains if the asym- 
metry summarized in Figure [14] is strictly a result of this 
one-sided asymmetry plus time-integration effects, or if an 
additional asymemtry is the cause of the blue and red knee in 
the fine-structure of Ha in LEO 16 and LEI 86, respectively. 

The "Bochum event" (Figure IT3T > has previously been in- 



terpreted as an asymmetric distribution of 56 Ni (e.g., iLucvl 
[19881: lHanuschik & Trumd[l990t IChugailll991ah . Although 
both blue and red emission features are observed in Fig- 
ure Q3J the blu e emission feature was only additionally ob- 
served in H/3 (lHanuschik & Thimml 119901) while the red- 
shifted em ission feature was a lso observed in infr ared hydro- 
gen lines dLarson et al.lll987t) and [Fe II] lines dHaas et alj 
119901) . IChugail d!991al) proposed a two-sided 56 Ni asym- 
metry: (1) a dominant 56 Ni cloud in the far hemisphere 
producing the red emission fine-structure and the redshifted 
emission peak at later epochs, and (2) a smaller, higher 
velocity 56 Ni cloud in the near hemisphere producing the 
blue fine-structure. However, instead of invoking a smaller 
56 Ni cloud in the near hemisphere, the more recent ex- 
planation for the blue fine-structure is a non-monotonic 
Sobolev optic al depth, r(v), for Ha with a minimum near v sa 
5000 km s" 1 dChugailll991tt lUtrobin et al.l ll995tlWanget"aI] 
I2002t lUtrobin & Chugail 120021 120051) . This interpretation 
does not require breaking sp herical symmetry to expla in the 
blue fine-structure. Although lUtrobin & Chugail (|2005j) favor 
this scenario, their model cannot successfully produce the op- 
tical depth minimum at the required strength. The smooth 
transition of the Ha profiles in PA from LE016 to LEI 17 
shows that the blue fine structure is dependent on viewing 
angle. As such, some form of deviation from spherical sym- 
metry must be the root of the blue fine-structure in the LE 
Ha profiles and presumably directly correlated to the blue 
fine-structure observed in the "Bochum event." 

The "Bochum event" was observed on days 20-100, while 
the redshift in lines of hydrogen and other elements was ob- 
served after ^150 days. For each LE spectrum, Table Q] lists 
the approximate range of epochs that contribute to the LE 
spectrum at greater than the 50% level. This relative con- 
tribution is based on the window function generated for each 
LE (Figure [3}, prior to any flux-weighted integration. While 
the temporal resolution in the LEs is not capable of distin- 
guishing between photospheric and nebular epochs, we can 
compare LEs with different degrees of nebular emission. It 
is therefore interesting to note that both the fine-structure and 
the velocity shift of the Ha line appear to be more dominantly 
a function of viewing angle rather than a function of epoch. 
LE066 probes the first ~ 80 days of the explosion only, but 
does not show fine-structure in Ha considerably different than 
the nearby LEs probing much later epochs. The fact that the 
blue fine-structure feature is also present in the early-epoch 
LE066 spectrum almost certainly links this feature to an ex- 
aggerated version of the original "Bochum event." 

Since the HaP-Cygni line is a blend with Ba II A64967, 
it is possible an asymmetry in the Ba II line is causing the 
observed blu e fine- structure in the LE spectra. However, 
lUtrobin et alj d 19951) determined that the inclusion of the Ba 
II line was not sufficient to explain the blue emission in the 
"Bochum event." Additionally, unlike Ha, the Ba II A6142 
line appears to be fit well with the isotropic model in the LE 
spectra. 

A two-sided 56 Ni model such as IChugail d!991al) could ex- 
plain the HaLE observations. A smaller high-velocity cloud 
blueshifted in the north causing the blue and red fine-structure 
in the north and south viewing angles, respectively. And a 
larger cloud redshifted in the south causing the red- and blue- 
shifted excess in the emission observed in the north and south 
viewing angles, respectively. The larger cloud could also 
be responsible for the Fe II asymmetry shown in Figure [T7] 
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The original observations of SN 1987 A provided strong evi- 
dence for the large cloud being in the far hemisphere as previ- 
ously stated. This is apparent in the isotropic model spectrum 
of LEI 86, where the time-integrated emission peak of Ha is 
shifted to the red. LEI 86 probes a larger range of epochs in 
the explosion (out to ~ 400 days), so the redshift and inte- 
grated line profile asymmetry is most apparent in the LEI 86 
isotropic model. The observed LE spectrum, however, shows 
a small blueshift with respect to zero velocity, implying the 
overabundance is inclined within 21° of the plane of the sky 
in the southern far hemisphere. If the cloud was ~ 20° from 
the plane of the sky, the radial velocities would be a factor of 
1.7—1.9 larger in the north and south viewing angles, more 
easily altering the Ha profile shape of the integrated LE spec- 
tra. 

In addition to SN 1987 A, a velocity shift in the Ha peak or 
a Bochum-like event has also been observed for the Type IIP 
SNe SN 1999em, SN 2000cb, SN 2 004di, and SN 2006oy , 
as well as the Type II SN 2006bc riElmhamdi et all 1 20031; 
Kleiser et al.ll201 It IChugai et al.ll2005t IChornock et al.lboiot 
Gallagher et all 120121). The Ha profiles of SN 2004dj are 
strikingly similar to that of the LEI 86 profile in Figure fT4l 
Although LE spectra cannot be compared dire ctly with in- 
stantan eous spectra, it appears the Ha lines from lChugai et al.l 
( 120051) would rese mble that of LEI 86 a fter being smoothed by 
time-integration. IChugai et al.l d2005l) were able to success- 
fully model the SN 2004dj Ha fine-structure using an asym- 
metric bipolar 56 Ni distribution, including a spherical com- 
ponent with two cylindrical components. They modelled the 
observer viewing a dominant jet 30° off-axis, resulting in a 
larger blueshifted emission peak with a smaller redshifted fea- 
ture in the profiles of Ha. This describes the LEI 86 LE 
profile in Figure [14] suggesting a bipolar 56 Ni distribution 
as a possible explanation. Since at the opposite viewing an- 
gle, LEO 16, we see the opposite set of features in Figure [Pfl 
the two LE viewing angles would have to be looking to- 
w ards opposite ends o f the bipolar distribution in the model 
of IChugai et all d200l . Since the viewing angles only differ 
by ~ 40°, this would put strict constraints on the orientation 
of the bipolar distribution. This is, however, consistent with 
the dominant 56 Ni component being aligned within ~ 20° to 
the plane of the sky, as previously suggested. 

As noted in the introduction, one of the benefits of LE 
spectroscopy is it allows the signatures of the explosion in 
the first few hundred days to be compared directly with the 
state of the remnant. The PA of the symmetry axis of the 
elongate d ejecta was meas ured to be 14° ± 5° using HST 
ima ging (IWang et al l 120021) and 15° ±0.9° using IFU spec- 
tra dKiaer et al.ll2010l) . This PA corresponds to LEO 16, where 
we see the maximum deviation from sym metry in Ha in the 
northern hemisphere. iKjaer et aTl (1201 Oh found the present- 
day ejecta to be blueshifted in the south and redshifted in the 
north, inclined out of the plane of the sky by ~ 25°. The 
bipolar 56 Ni distribution proposed above to explain the LE 
observations is therefore aligned with the roughly 25 year-old 
ejecta both in PA and inclination out of the sky. The inner 
ring of circumstellar material is i nclined 49° out of the plane 
of the sky, redshifted in the north (ISugerman et al.ll2005l) . The 
outer circumstellar rings are similarly inclined, presumably 
related to the rotation axis of the progenitor. The elongated 
ejecta and 56 Ni distribution probed by the LE observations are 
therefore aligned approximately in plane with the equatorial 
ring, as opposed to sharing a symmetry axis as initially pro- 



posed by IWang et al.l ( 120021) . This disfavors the axially sym- 
m etric jet-induced ex plosion model for SN 1987 A proposed 
bv lWang et all ( 120021) . 

The "mystery spot" of SN 1987A was a bright source ob- 
served in speckle interferometry measurements 30, 38, and 50 
day s after the explosion, separated by only 6 mas from the 
SN (iNisenson et a l.| [l987t iMeikle et al.|[l987t). and appearing 
at PA 16°/ 196° ± 2°. INisenson & Papalioliosl (119991) repro- 
cessed this data, and identified two mystery spots in the data 
at opposite position angles, along the same 16°/ 196° axis. 
In order to be associated with the SN, they claimed the mys- 
tery spots must be at relativistic speeds with the northern spot 
blueshifted. A satisfactory explanation of the mystery spot (or 
spots) has yet to appear, a fact which is often forgotten. Since 
the PAs for the LE Ha spectra with maximum fine-structure 
match the mystery spot PAs to within 10°, future modelling 
of the LE lines may aid in an explanation. 

If in fact the blue fine-structure of the "Bochum event" is 
due to an asymmetrical 56 Ni feature in the northern hemi- 
sphere, as the LE spectra here suggest, 56 Ni is mixed to even 
higher velocities than previously considered. The blue feature 
emerges after 20 days at a radial velocity of -5000 km s" 1 in 
the Ha profile. Since the blue fine-structure is most domi- 
nant in LE016 at a scattering angle of 17°, the lower limit 
of the absolute velocity of the 56 Ni cloud in the near hemi- 
sphere is 5200 km s" 1 , a velocity difficult to explain in cur- 
rent neutrino-driven explosion models. Absolute velocities 
of > 7000 km s" 1 are required if the cloud is inclined within 
45° to the plane of the sky and upwards of 10000 km s" 1 
if within 30°. Assuming the dominant southern 56 Ni over- 
abundance (and change in Ha peak velocity) is correlated to 
the red emission feature observed in the "Bochum event" re- 
quires absolute velocities of ~ 10000 km s" 1 for that feature 
for an inclination of 20° out of the plane of the sky. For the 
case of SN 2000cb. lUtrobin & ChugaHd201 ll) required mixing 
56 Ni to velocities of 8400 km s" 1 to reproduce the observed 
lightcurve of SN 2000cb, although it was a more energetic 
Type IIP than SN 1987 A. Therefore, these unrealistically high 
velocities of 56 Ni for current neutrino-powered core-collapse 
explosion models require us to step back and view the asym- 
metries in Figure[l4]at their most basic leve l: a non-spherical 
excitat ion structure in the early ejecta. [Dessart & Hillierl 
d201 ll) demonstrated that the P-Cygni line profiles from non- 
spherical Type II ejecta can be altered significantly depending 
on viewing angle. Only future modelling of the LE spectra, 
taking into account the effects of time-integration, will allow 
the ejecta geometry to be determined or further constrained. 

5. CONCLUSIONS 

We have obtained optical spectra from the LE system of 
SN 1987 A, allowing time-integrated spectra of the first few 
hundred days of the original explosion to be viewed from mul- 
tiple lines of sight. We have modelled the LE spectra using 
the original photometry and spe ctroscopy of SN 1987 A using 
the model of iRest et al.l (1201 lal) . Using specific examples we 
have showed the model correctly interprets LE flux profiles 
and spectra when both scattering dust properties and observa- 
tional properties are taken into consideration. 

Each PA on the LE system represents a unique viewing 
angle onto the photosphere which can be spectroscopically 
compared to the isotropic spectrum calculated from the origi- 
nal outburst of SN 1987 A. The LE spectra show evidence for 
asymmetry, demonstrating the technique of targeted LE spec- 
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troscopy as a useful probe for observing SN asymmetries. The 
observed asymmetries can be summarized as follows: 

1 . Fine-structure in the Ha line stronger than the original 
"Bochum event" is observed most strongly at PAs 16° 
and 186°, with the Ha profiles showing opposite asym- 
metry features in the north and south viewing angles. 

2. At PA 16° we observed an excess in redshifted 
Ha emission and a blueshifted knee. At PA 186° we 
observed an excess of blueshifted Ha emission and a 
red knee. This fine-structure diminishes slowly as the 
PA increases from 16° to 16° + 90°, with the LEs per- 
pendicular to the symmetry axis showing no evidence 
for asymmetry. 

3. At PA 16° we observe a velocity shift in the Fe II A5018 
line of the same magnitude and direction as the veloc- 
ity shift in the Ha emission peak at the same viewing 
angle. As with the Ha fine-structure, this velocity shift 
appears to diminish as PA increases from 16°. 

4. At PA ~ 1 15° (roughly perpendicular to symmetry axis 
defined by 16°/186° viewing angles) we observed an 
unidentified feature near A5265A not observed in any 
other viewing angle. 
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This symmetry axis defined by the 16°/ 186° viewing an- 
gles is in excellent agreement with the current axis of symme- 
try in the ejecta geometry, the initial polarization and speckle 
observations, as well as the location of the "mystery spot." 
The Ha lines at PAs 16° and 186° are very similar to the 
Ha lines observed in SN 2004dj a n d mod elled as a two-sided 
56 Ni distribution by IChugai et alJ (120051) . This same model 
could describe a two-sided ejection of 56 Ni in SN 1987 A as 
probed by the LE Ha lines. The 56 Ni is blueshifted in the 
north and redshifted in the south, with the dominant overabun- 
dance of 56 Ni being inclined ~ 20° from the plane of the sky. 
While these observations strongly argue for a two-sided dis- 
tribution of high-velocity 56 Ni, only future modelling of the 
Ha line-formation will be able to confirm this. The indication 
that high-velocity 56 Ni is not correlated with the inner ring 
and presumed rotation axis may indicate that the explosion 
mechanism is independent of rotation. 
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APPENDIX 

EFFECT OF DUST SUBSTRUCTURE IN LE PROFILES 



For the observed LE profiles presented in this work, 11/14 show no strong evidence for dust substructure within the profile. 
These LE profiles can be modelled with one dust sheet effectively. However, LEI 13, LE180, and LE186 all show significant 
substructure in the their profiles (Figuresl4"land[To*l) and we seek to quantify the effect this can have on the observed LE spectrum. 
Figures [T9ll2"Tl show the observed LE profiles and the effect of dust substructure on the window functions and effective lightcurves. 
Figure [22] shows the differences in Ha strengths obtained with and without including the effect of the secondary dust structure. 
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FIG. 19. — LEFT: LE113 observed LE profile (black points) with the best fit to the primary dust sheet (solid red line), as well as the best fit to the secondary 
dust sheet (dashed red line). The location of the slit is shown with the grey region. MIDDLE: Window functions for the primary (solid red line) and secondary 
(dashed red line) dust sheets, as well as the total window function (solid black line). RIGHT: Effective lightcurves for the primary dust sheet (solid red line), the 
total effective lightcurve from both dust sheets (solid black line), and the original lightcurve of SN 1987 A (dotted line). The effect of the secondary dust sheet is 
minimal in this case. 
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FIG. 20. — LEFT: LE186 observed LE profile (black points) with the best fit to the primary dust sheet (solid red line), as well as the best fit to the secondary 
dust sheet (dashed red line). The location of the slit is shown with the grey region. MIDDLE: Window functions for the primary (solid red line) and secondary 
(dashed red line) dust sheets, as well as the total window function (solid black line). RIGHT: Effective lightcurves for the primary dust sheet (solid red line), 
the total effective lightcurve from both dust sheets (solid black line), and the original lightcurve of SN 1987A (dotted line). Without the secondary peak, the 
effective lightcurve is significantly biased towards later times due to the location of the slit. The secondary peak contributes flux from early epochs, reducing the 
late-epoch biased produced by the slit. 



As shown in Figure [22] the effect of the secondary substructure on the resulting Haprofile is minimal for LEI 13 (< 4% 
difference in strength of emission peak), primarily because the flux contribution from the wing of the LE profile is so small. 
For LEI 86, including the secondary LE profile decreases the Ha emission peak by 11%. Note, however, that because of the 
significant offset of the slit from the primary peak, the Ha emission for LEI 86 is in excess of the full lightcurve- weighted 
integration of SN 1987 A. For the case of LEI 80, the secondary peak significantly alters the window function and subsequent 
Ha line, increasing the emission strength by ~ 40%. 
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FIG. 21. — LEFT: LE180 observed LE profile (black points) with the best fit to the primary dust sheet (solid red line), as well as the best fit to the secondary 
dust sheet (dashed red line). The location of the slit is shown with the grey region. MIDDLE: Window functions for the primary (solid red line) and secondary 
(dashed red line) dust sheets, as well as the total window function (solid black line). RIGHT: Effective lightcurves for the primary dust sheet (solid red line), 
the total effective lightcurve from both dust sheets (solid black line), and the original lightcurve of SN 1987A (dotted line). Here the two dust sheets create LE 
substructure with nearly equal amounts of flux. The window function in this case is significantly shifted to later times due to the presence of the secondary peak 
at larger values of p. 
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FIG. 22. — Resulting Ha profiles obtained by integrating the historical spectra of SN 1987A with the effective lightcurves in Figures IT 912 1 1 for LEI 13 (left 
panel), LE186 (middle panel), and LE180 (right panel). For each panel, the solid and dashed lines are from integrations of the total and primary-only effective 
lightcurves. The dotted line is Ha profile obtained by integrating the actual lightcurve of SN 1987A with no window function truncating the integration. 

THE EXCESS OF Ha EMISSION IN LE186 

The observed LE spectra and isotropic models for LE180 and LE186 in Figure QT| show two different results. LE186 shows 
a clear, high signal-to-noise excess in Ha emission over the isotropic dust-modelled spectrum. LEI 80, probing essentially the 
same viewing angle onto the photosphere, does not show such an excess in the Ha profile. In addition to the much lower signal- 
to-noise of LEI 80, we can use arguments based on the above dust substructure analysis to show that the excess emission observed 
in LEI 86 is physical. 

As shown in the middle panel of Figure [22] the effect of the secondary LE peak in the LEI 86 profile is to reduce the emission 
peak of Ha . The maximum amount of Ha emission for LEI 86 therefore occurs when no secondary substructure is included. 
Since this gives an excess of only ~ 10%, even this scenario is not able to account for the ~ 30% excess in Ha in the observed 
LEI 86 spectrum. 

For the case of LE180, the effect of the substructure is to increase the Ha emission peak substantially (~ 40% using the profile 
fits in Figure |2iTi. Unlike the case of LEI 86, where the LE profile fits are well constrained within the data, the fits to the LEI 80 
profile in Figure[2T]are an example of a situation where our profile fitting algorithms cannot extract meaningful results due to the 
intrinsic complexity and low signal-to-noise ratio. It is most likely the case the secondary substructure in the LEI 80 profile is 
actually much wider than shown in Figure |2T] which would result in more late-epoch contribution from the secondary peak and 
more emission in Ha . 



